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Abstract 

We calculate the contribution of the weak annihilation to the B — » /ry 
decay amplitude by means of QCD sum rules using the photon light-cone wave 
function. We find that this long-distance contribution amounts to about 10% 
of the leading short-distance effect in B~ — > p~7- On the other hand, weak 
annihilation is the dominant source of the corresponding D meson decays and 
according to our estimates, yields branching ratios of O(10 -4 ) for D° — ► K*°j, 
O(10~ 5 ) for D s -» p+7, and O(10~ 6 ) for D~ -> p~7 and for D° -> p°7- 
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1. Rare decays of 5-mesons, such as the recently observed processes B — > K*<-y 
and B — > X s 7 j[J |2| are becoming an important tool for studying new forces beyond 
the standard model ||. Indeed, many authors have investigated the effects of multi- 
Higgs models, supersymmetric theories, left-right symmetric models, etc. on this 
process The interest in these decays stems from the fact that they occur only 
through loops and are therefore particularly sensitive to "new physics" . 

Besides the b — > s transition, there are also b — > d processes which are suppressed 
by the ratio V t d/V ts ~ 0.1 of CKM matrix elements. Within the standard model 
the rates are then estimated to be about 10~ 5 and 10 -6 for b — > d^y and B — > fry, 
respectively. This mode is sensitive to other CKM matrix elements (and possibly to 
other new physics) and may show large CP- violation ||. 

The "inclusive" processes like b — > 57 can be usually evaluated on the quark 
level within perturbative methods [], but are more difficult experimentally. Exclusive 
processes, on the other hand are more easily observable, but require difficult non- 
perturbative calculations of matrix elements in order to yield useful results. 

It is generally believed that the short distance penguin mechanism dominates 
the exclusive decays. The corresponding matrix element has been calculated with 
various methods |5], |6| which yield similar results; this indicates that the related 
uncertainties are presumably under control. Besides the matrix element of the pen- 
guin operator, there are other long-distance effects which must be analyzed. Several 
authors have recently investigated these effects 0, §, ||, |10|, [H], Q in a phenomeno- 



logical way. The crude estimates show that the long distance contributions can 
indeed be small. However, the results obtained vary enormously and can give at 
best an order of magnitude estimate. Clearly, a more reliable calculation is desired. 

We may distinguish two long distance effects. First, there is a penguin type 
mechanism where the low momentum part of the usual penguin loop must be in- 
cluded. The above mentioned estimates so far indicate that it is at most several 
percent of the short distance result both for exclusive and inclusive processes. In 
these estimates, the radiative transition B —>■ M7 is modelled by a sum over the 
processes B — > MV* followed by the transition V* — ► 7 where V* is a virtual vec- 
tor meson with appropriate quantum numbers. The problems with this picture are 
obvious: the effective couplings (for the primary vertex and for the transition of 
the resonance into a photon) are only measured at the mass of the resonance and 
must be scaled to zero. Then, some information about the various structures of 
the primary vertex must exist, since only the transversal part can contribute to the 
photon transition. In the existing literature, one or the other of these problems is 
not addressed, but there is generally a strong suppression. For the particular case 
of B — > if* 7 this suppression however, crucially depends on the model assumptions 
about the share of the factorizable contribution. If factorization dominates, we can 



combine the results of refs. |12| and |J to predict a negligible effect for the J/^f 
resonance. But only a more systematic calculation combined with direct estimate 
of nonfactorizable effects can clarify the role of this mechanism. 



1 An exception is the "long distance" part of the penguin mechanism (see below) 
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Then, there is a completely different class of diagrams which do not involve 
loops of heavy (new) quarks but just the ordinary four fermion interaction. In these 
graphs, the 'spectator' participates in the weak annihilation; we therefore call them 
weak annihilation (WA). Because of the CKM matrix elements, this mechanism is 
negligible for B — > K*j but important for other exclusive radiative decays. For 
definiteness, we will consider it for the decay B — > pj, although other decays can be 
treated in similar ways and will be discussed briefly at the end. Fig. 1 represents 
schematically the mechanism. A perturbative evaluation of the diagrams in Fig.l 
0, |10| is rather uncertain, due to the almost on-shell propagation of the light quark 
which implies the use of a the poorly understood 'constituent' quark mass. The 
phenomenological description via the chain B — > pp* — > pq |I| suffers from the 



unknown Bpp coupling and the problems connected to the off-shellness of the p. 

In this paper we show that existing QCD sum-rule technology allows us to provide 
reliable predictions for the long distance effects. We will discuss in this short note 
only the WA effects which are more uncertain in the phenomenological approach and 
which are important by themselves to a large class of transitions. We will only be 
interested in the leading effects; more accurate and further results will be presented 
in a future publication; we will asses here some of the remaining uncertainties. 



2. The relevant effective Hamiltonian for B — > /ry consists of two operators [p~3 1 : 



H w = -^=V ub V: d {( Cl (dL, lU )(uL%) + c 2 {uL^u)(dm)} , (1) 

where = 7^(1 — 75), and c\ and c 2 are the Wilson coefficients. For the decay 
B~ — > p~7 the combination 

H w = ^Vu^MdL^iuL^b) , (2) 

with a\ = Ci + C2/3 enters. Similarly, the combination a 2 = Ci/3 + c 2 multiplies 
the corresponding operator for B° — > p 0/ ~f. Following the phenomenological ansatz 
|L4|, the coefficients a±, a 2 are extracted from two-body non-leptonic decays using 



factorization. Recent fits |L5[ indicate a considerable departure of a 2 from its short 
distance value (see, however, the recent discussion in The reason for that 

departure may be connected with nonfactorizable contributions to the hadronic ma- 
trix element, which in the decay B~ — ► p~7 come from the color octet operator 
(dL^u)(uU L ^b) with coefficient 2c 2 . 

The actual size of nonfactorizable effects is under discussion at present. First 
direct estimates [|17]], |18|,p9| by QCD sum rule methods indicate that nonfactoriz- 



able effects in B decays are quite process dependent. Moreover, it is not clear if one 
can use the same value of the effective parameters a 12 for nonleptonic and radiative 
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decays, processes with different physical final states. Nevertheless, since the decay 
B~ — > depends on the parameter a\ which is near its perturbative value a\ <~ 1, 
we assume for the present work that factorization is sufficient. At the same time 
one should keep in mind that the result for B° — > p°7 which depends on 0,2 will only 
be approximate. 

Turning to more detailed analysis of the WA mechanism we first of all neglect 
the photon emission from the final state quarks invoking the well known helicity 
arguments: being factorized, the matrix element of (f2|) is in this case proportional 
to the light quark masses: (p~ \ Hw \ -£>~) 7 ~ P^bIb{p~ | (dL^u) | 0) 7 ~ 0{m Ul md). 
Note however that non-factorizable interactions such as gluon exchanges between 
final and initial quarks lift this suppression. 

Thus we are left with two diagrams in Figs. la,b where the photon is emitted 
from the initial b and u quark lines. The corresponding matrix element can be 
written as 

(p~ I H w I £T> 7 = ^=V ub V: dai f p m p ^ I (uL,b) | B~) 7 , (3) 

where we used (p~ \ (dL^u) | 0) = f p m p e^, denoting by f p and e p the decay constant 
and the polarization vector of the charged p-meson. It then remains to calculate the 
matrix element 

(0 I (uWb) I £-} 7 = -A PC e prX(7 p T e x q r7 + iA PV [q^e ■ p) - e p {p ■ q)] , (4) 

which describes the annihilation of B~ into the current uL^b with momentum p 
accompanied by the emission of a real photon with momentum q and polarization 
vector e. We have written explicitly decomposition of this matrix element into 
two gauge invariant (with respect to the electromagnetic field) terms Apc(Apy) 
corresponding to parity conserving ( parity-violating) B —>■ p7 decay amplitude. 
We note that for the short distance penguin mechanism one has Ape = Apy if the 
light quark mass is neglected ( see e.g. [[J). In principle, polarization experiments 
could distinguish the two amplitudes and thus determine to what extent this equality 
is valid. In fact, it turned out that the WA mechanism considered here does not 
respect it, but the deviation is due to nonleading contributions and consequently 
small. 

As noted, we propose to use QCD sum rules in order to calculate the matrix 
element (0). Since the photon emission from the light quark takes place at large 



distances, the use of standard QCD sum rules based on the local operator 
product (OPE) expansion is not sufficient. Rather, one should use a light-cone 
expansion which is adequate for exclusive processes with light particles. It will 
involve the hadronic wave functions on the light-cone which encode the photon 
emission by a quark- ant iquark pair at light-like separation in close analogy to the well 



known pion wave function The photon light-cone wave function was introduced 



in ref. [22j for calculating the amplitude of weak radiative decay S — > p7 and used 
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later in ref. P3) to evaluate the nucleoli magnetic moments. 



3. We now proceed to calculate the correlation function 

n„(p,g) = i J d 4 xe l P x (0\T{u(x)L^b(x),m^u(0)}\0)F( q ) ■ (5) 
in the external electromagnetic field 

F a p{q, x) = i(epq a - e a qp)e iqx (6) 

with momentum q and polarization vector e. The function 11^ can be decomposed 
into two parts corresponding to parity conserving and a parity violating part in (|4]): 

n M = Upc4e flTXp p T e x q p + Tl PV [q^e ■ p) - e^p ■ q)] (7) 

In the region (p + q) 2 < and at p 2 = m 2 << m 2 , the heavy 6-quark is far off-shell. 
In particular, photon emission from the heavy 6-quark takes place perturbatively. 
The accompanying light-quark propagator may then be described by the local OPE 
containing as a first approximation the free propagation and in the next orders, 
interaction with quark and quark-gluon vacuum condensates. The corresponding 
contributions to the correlation function are depicted in Figs. 2a, b and c, re- 
spectively. As far as photon emission from the light w-quark is concerned, after 
contracting the 6-quark line, one is left with matrix element 

n„(p,9) = i J (2 ^ ) f ( ^ fc2) e- (p - fc) -(0|^(x)L M K+ ^) 75M (0)|0) F(9) . (8) 

The diagram Fig. 2d describes only the short-distance part of this matrix element 
corresponding to the photon emission from freely propagating u quark. 

To take into account the long-distance part one needs to introduce additional 
nonperturbative parameters, describing the light-quark propagation in the presence 
of external electromagnetic field which is schematically shown in diagram of Fig. 
2e. In the case of constant electromagnetic field (q — > ) one still can use the 



local OPE i.e. expand (gj) around x = in the external field [24, 25]. Physically, 
the most important parameter emerging in this expansion is the so called induced 
quark condensate 

{0\qa a/3 q\0) F = e q (qq)xF af3 (9) 

where e q is the quark electric charge, (qq) is the quark condensate density and pa- 
rameter x has physical meaning of magnetic susceptibility of the quark condensate. 

However, here we deal with essentially large photon momenta of the order of 
heavy quark mass. In this kinematical configuration the use of local OPE will 
lead to unmanageable infinite series of condensates for the expansion of (B) (see a 
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recent detailed discussion in P5| ). The adequate approach is to use the expansion 
near the light-cone x 2 = introducing matrix elements of nonlocal operators. The 



leading twist two contribution |22| , |23| emerges from the expansion of the operator 
u(x)a a pu(Q) and the corresponding part of the matrix element in the external photon 
field can be parametrized asf] 

(0\u(x)cr a(3 u(0)\0) F{q) = e u (uu) f du(p(u)F aP (ux) . (10) 

<J 

Here, the function (p(u) has the meaning of the photon wave function in terms of 
its quark-antiquark constituents and may be interpreted as the distribution in the 
fraction of light-cone momentum q + q 3 of the photon carried by a quark. The 
asymptotic form of this wave function is known |2l], |22j: 

ip(u) = 6xu(l - u) , (11) 

where the appearance of the magnetic susceptibility x is evident from taking local 
limit of ( |Tt]| ) and comparing with the definition ([|). Moreover, according to the 
analysis carried out in ref. p2| , non- asymptotic effects in <p(u) and higher twist 
(twist 4) contributions to (|T0|) are small, contrary to the case of the pion wave 
functions. Therefore, in what follows we will use the asymptotic wave function (|TT|). 
A rough estimate of the twist four contribution carried out below will allow us to 
neglect all higher twist effects. 

To proceed, we retain in eq. (||) the parts containing u(x)a tlu u(0), insert eq. 
([T0|), and after integration over x and k obtain the leading twist 2 contribution to 
the correlation function ^) corresponding to the long distance photon emission from 
u-quark and depicted in Fig. 2e: 

ngua = n « = l6u{Uu) f 1 du (12) 

Jo mi — (p + uq) z 

The next important contribution stems from the perturbative loop diagrams of 
Figs. 2a and 2d. Their contribution to the parity conserving invariant amplitude 
written in the form of dispersion integral in variable (p + q) 2 at fixed p 2 is: 



t 3im h r°° ds 



'-PC 



An 2 J ml [s — (p + q) 2 ](s - p 2 



. s - ml s 

[e u - e b ) h e b ln— 5 

s mi 



(13) 



where the pieces proportional to e b and e u correspond to the diagrams Fig. 2a 
and 2d, respectively. The contribution of the loop diagrams to the parity violating 
amplitude is: 



yrpert 3i?7lh f°° ds 
11 1 



L PV 



An 2 J m \ [s — (p + q) 2 ](s — p 



2\2 



\m 2 h -s-t< 



2 the path-ordered exponential gauge factors for both gluon and photon fields are not shown 
for brevity 
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e b I 2m 2 b - s - — ^ 1 }- ^ + e b (2m z b - s +p 2 )ln-^ 



m' b 



(14) 



We see that the amplitudes Upc and Upv are indeed different. To assess the impor- 
tance of nonperturbative effects in the case when photon is emitted from the heavy 
quark we have also calculated the contributions of quark condensate corresponding 
to Fig. 2b with the result 

n m _ n <??> _ W 

( m b -p )i m b - (p + q) ) 

In calculating the perturbative and quark condensate contributions, the contact 
terms violating gauge invariance with respect to the photon field should be properly 
taken into account. If one carries out the calculation in the fixed point gauge for 
the photon field, the origin must be chosen carefully. The quark-gluon condensate 
contribution (Fig. 2c) and higher-twist contributions corresponding to three-particle 
photon wave function (Fig. 2f) are neglected in the present calculation. 

To understand the level of accuracy for the most important contribution con- 
taining the photon wave function we need some estimate of nonleading higher twist 
contributions neglected in the matrix element in eq. (0). We parametrize the 
next-to-leading twist 4 contribution in the expansion of this matrix element as : 



\u(x)a al3 u(0)\0) t F v ^ ) u = e u (uu) / du{/i(w)F a/3 (ux)a; 2 

J 

+f 2 (u) [F ap (ux)x p x l3 - F ljp (ux)x p x a }} . (16) 



As a result, 



Uf' sU = -4ie u (uu) [ du- 
Jo 



fi(u) 



mt 



2m 2 h 



— (p + uq) 2 } 2 ml — (p + uq) 2 



(17) 



and the corresponding answer for n^y sM is obtained by replacing fi(u) in the numer- 
ator of fllTD by the combination fx + f 2 . The only information about distributions 
/i 2 is available from comparison with the twist 4 terms in the light-cone sum rules 
for the S — > jry decay p2] and for the nucleon magnetic moment |23| which result 
in the following relation: 

fi(u) + ^ = -\[(l-u)-Mu)} ■ (18) 

The first term in r.h.s. of ([Tj|) originates from the expansion of the light quark field 
in the electromagnetic field, whereas the wave function <^ 4 corresponds to analogous 
expansion in the gluon field and is related to the three-particle photon wave functions 
of twist four via QCD equation of motion. According to estimates of [^|, <p4,(u) ~ 
30ku 2 (l — u) 2 + with the parameter k ~ —0.2 and the nonasymptotic corrections 
denoted by ellipses and not shown for brevity. We assume for an order of magnitude 
estimate that fi = f 2 , having in mind that the normalization integrals of these two 
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distributions are equal. We plan to analyze this contribution in more details in 
future. 

4. The QCD sum rule is obtained as usual by equating the hadronic representa- 
tion of the correlation function 11^ with the result of the QCD calculation, a sum of 
all contributions calculated above. Inserting hadronic states with B meson quantum 
numbers into eq.([|), we have the following decomposition for invariant amplitudes 
in eq.(0): 

ifBm 2 B A PC (pv) 



n 



00 ds P h pc(PV)( s >P 2] 



m b [m B 2 - (p + q) 2 \ J so s - (p + q) 2 

where the first term represent the I?-meson contribution containing the matrix el- 
ement (H) and (B \ bij 5 u | 0) = fBm 2 B /mb. The second term in (|l^) represents 
contribution of the higher states in B channel starting from some effective thresh- 
old so- We invoke the usual hadron-quark duality and replace the spectral density 
Ppc(pv) by the imaginary part of Upc(pv) calculated in QCD. Equating now eq. 
( |T9| ) with the result of this calculation given in eqs. (|T^)-(|T5|) and applying a Borel 
transformation in (p + q) 2 to suppress the higher states we obtain the desired QCD 
sum rule for both amplitudes entering the hadronic matrix element (^). The result 
for the parity conserving amplitude is 



.4 



PC 




1 du 



fB.m 2 



u 



exp 



M 2 



p 2 (l 



u 



uM 2 



e u {uu)(p(u) 



((«._ 

An 2 \ mi 



p 2 )(l 



p 2 (l — 11 



e b (uu) 



mi 



p A 



exp 




p\\ 



umi 



and for the parity violating amplitude 
mj I f l du 



A 



PV 



fprn 2 




u 



exp 



ni 



B 

M 2 



mi 



p 2 (l 



u 



uM 2 



u 



e u {uu)ip(u) 



(20) 



+ e b ln 




p 2 u 2 



m b — P 2 (l — m)/ 

{ml — p 2 )(l — u) 
ml — p 2 {l — u) 



2 2 ex P 
mi — p z 



1 1 i s 
m B — mi 



M 2 



(21) 



where M 2 is the Borel parameter. 

In order to compactify these expressions, the dispersion integrals for the loop con- 



tributions are written in terms of a new integration variable u 



[mi 



■p 2 )/(s 



P 2 ), 



so that substitution of A 



(mi 



p )/(sq — p ) instead of the lower limit in all 
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integrals in (|20D,(|2TD is equivalent to subtraction of higher state contributions. Note 
that the results of our calculation are, in general, the two form factors Apc,pv(p 2 ) 
determining p 2 dependence of the matrix element (f|) . For the problem under inves- 
tigation we need only the values of these form factors at p 2 = m 2 or at p 2 = m 2 K ». 
The sum rules are valid in a much wider region of timelike variable p 2 , paramet- 
rically up to 0(m 2 — 0(GeV 2 )) and practically up to p 2 =15 GeV 2 (see analysis of 



analogous sum rules for heavy-to-light form factors in B 26, 27 1). 



5. We proceed now with the numerical analysis. The important point is that 
all parameters entering the l.h.s of eqs. (p0|), (pl|) are known since they also enter 
various other QCD sum rules; the universality of nonperturbative inputs is of course 
the main advantage of this approach. The value of magnetic susceptibility was 
determined several times [ 28[ , ||29|| with essentially the same result x = — 4.4GeV -2 



at the normalization scale of 1 GeV. We use (uu) = —(240 MeV) 3 at this scale. 
The anomalous dimension of the qo^q operator is known |30| to be (-4/27) whereas 
the anomalous dimension of the quark condensate operator is (4/9) resulting in 



X = —3.4 GeV~ 2 and (mm) = —(260 MeV) 3 at the scale p^ ~ yfnp ~ m l inherent 
to the correlation function (^|). The comparatively large numerical value of the 
parameter x justifies the use of leading twist approximation. The relatively high 
scale Hb is an additional argument on favor of the asymptotic form ( |TTD throughout 
our numerical calculation. 

The values of all parameters corresponding to the 5-channel are known from 
two-point sum rules for the b-quark current: fp = 140 MeV, m h = 4.7 GeV, s = 
35 GeV 2 . (see e.g. p6j, |27|). Note that the £?-meson decay constant is taken without 



0(a s ) corrections since they are also not included in the sum rules (^) and ( pTj ) 
where they partially cancel those of fp. With these parameters, the amplitudes Ape, 
Apv are calculated at p 2 = m 2 as a function of the Borel parameter. The predictions 
of the sum rules are very stable in a rather wide region of Borel parameter and vary 
by at most 5% with the changes of nib, ^o, fs within the intervals allowed by the 
two-point sum rule for fp- In analogy with other sum rules we isolate a certain 
interval of Borel parameter which we determine as 8 < M 2 < 12GeV 2 . Here our 
estimate of twist four operators gives contributions below the 5% level and higher 
states according to the quark-hadron duality model contribute less than 30%. Our 
final prediction for the amplitudes defining the matrix element (|4]) for the WA decay 
B~ — > p~7 is : 

A PC (B- -> p-j) = 1.1 • lO^GeV" 1 , A PV (B- -> p~j) = 0.85 • lO^GeV" 1 . (22) 

with negligible variation within the interval of M 2 specified above. The correspond- 
ing prediction for the amplitude of the neutral mode B° —>■ p°7 is easily obtained 



S 



by replacing u — > d in the sum rules. The prediction in the same Borel window is : 

A PC (B° -> p° 7 ) = -7 • Kr 3 GeV _1 , A PV (B° -> p° 7 ) = -5 • lO^GeV- 1 . (23) 

All these predictions have a negligible variation with the Borel parameter within the 
interval specified above. As expected, an important part of the amplitudes (^2|) and 
(f23|), (40^-50)% comes from the contribution of the photon twist-two wave function. 
The contribution from the loop diagram with photon emission from the light quark is 
about the same and amounts to (50-r-60)% of the total result. The photon emission 
from the heavy quark (both perturbative and nonperturbative) is negligibly small. 
The strong isospin violation (or in other words, dependence on the flavour of the 
spectator quark) which is evident from the ratios of the amplitudes of charged (52]) 
and neutral ( ^3"D modes, is one of the characteristic features of the WA mechanism 
in contrast to the short distance penguin amplitudes which are independent of the 
flavour of the spectator quark. 

Concerning the accuracy of our calculation we first notice that omitted higher 
order terms are indeed very small. The higher order quark-gluon condensate con- 
tribution, Fig. 2c, can be safely neglected since the quark condensate contribution 
itself turns out to be at the level of 1%. Also the estimate of twist 4 contribution al- 
lows to neglect the diagram Fig. 2f containing the twist 4 three-particle photon wave 
function. One may even argue that this diagram is additionally suppressed because 
the gluon is emitted from the heavy quark line (such a suppression is present in other 
light-cone sum rules for the heavy-to light transition amplitudes, see e.g. [^7|, |2(J ). 
The higher twist terms and nonasymptotic corrections are the main source of uncer- 
tainty in the light-cone sum rules involving the pion wave functions. In the present 
situation, the main uncertainty comes from the range of the parameter x which is 
known with 10 % accuracy. Our conservative estimate of the overall accuracy is 
about 15-20 %. Furthermore, there are 0(a s ) corrections to the correlation function 
d^). Presumably, the largest part originates from perturbative gluon exchanges in 
the leading diagram of Fig. 2e. Their calculation is straightforward, but beyond the 
goal of our work. 

It is immediate to convert eqs. (^) and ( |2"T| ) into the sum rules for the analogous 



hadronic matrix elements (0 | (d(u)L^c) | _D+^) 7 determining the radiative decays 
of charmed mesons D —> fry via WA mechanism. One has just to replace b with 
c , B with D in the sum rules (|20|), ( pT|) and substitute corresponding parameters: 
f D = 170 MeV, m c = 1.3 GeV, s = 6 GeV 2 , (we take them from fH) X = -4-0 
GeV" 2 and (uu) = -(250 MeV) 3 . 

In the Borel parameter range 3.0 < M 2 < 4.5 GeV 2 which satisfy all usual 
requirements we obtain (with small variations) 

A PC {D + -> p+ 7 ) = -1.7-10- 2 GeV-\ A PV (D+ -> p+ 7 ) = -1.5-10- 2 GeV -1 . (24) 

A PC {D° -> p° 7 ) = +9.6 • 10- 2 GeV-\ A PV (D° -> p° 7 ) = +5.5 • lO^GeV" 1 . (25) 

For the correlation function (f|) with c quark currents the hierarchy of contributions 
drastically changes. The magnitude of photon emission from the c quark is almost 
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as high as from the light quark. At the same time the relative role of the nonpertur- 
bative emission from the light quark (i.e. of the photon wave function contribution) 
compared with the perturbative emission is larger. As a result, the amplitudes ( f24|) 
and ([25|) represent interferences of various important contributions and the role of 
subleading terms is in general higher, especially for D~~ decay where the estimate of 
the size of the twist 4 term indicates a 20% correction. For this particular mode a 
higher accuracy is required, including the calculation of terms omitted in our present 
calculation. 

For completeness we list the corresponding amplitudes for the remaining B and 
D modes which are simply connected to the modes considered above : 

A(B~ -> K*~i) ~ A(B~ -> p~7) (26) 

and 

A(D° -> ir° 7 ) ~ A(D° -> p°j), A{D + -> p + 7 ) ~ A(D S -> p+j). (27) 

These relations are valid for both PC and PV amplitudes in the SU(3) flavour limit, 
i.e. if one can neglect the variation of form factors Apc,pc from p 2 = m 2 to p 2 = m 2 K * 
or, in the case of the second relation in (p7|), the difference between photon emission 
from s and d quarks. From our sum rules we are able to safely predict that the 
amplitudes at the point p 2 = m 2 K , are 0(1%) ( 0(10%) ) higher for B ( D). The 
WA mechanism is absent for the modes B° —>■ K* 0/ y, B s — > K*°j,B s —>■ 07. and for 
D~ — ► K*~j) (with the possible exception of CKM suppressed contributions). 



6. With the amplitudes determining the heavy meson annihilation matrix ele- 
ment we can finally compare the WA effects to that of the short-distance penguin 
mechanism and obtain the corrections to the branching ratios. 

The branching ratio corresponding to the decay amplitude (0) with the photon 
transition matrix element defined in is: 

/ 2 2 \ ^ 

BR(B- - p-i) = ^{V ub V: dai f f p m 2 p [ ^" j [Ale + A 2 PV ]r B - (28) 

for B~ — > p _ 7 and with obvious substitutions for the other decays. Taking for the 
present analysis ai = 1 , f p — 216MeV from the leptonic width of p-meson which 



we take along with all other parameters entering (|28| ) from JHJ] we finally obtain for 
the branching ratio corresponding to the WA mechanism : 

BR(B- - p- 7 ) ~ 7.0 ■ 10- 5 GeV 2 (^) ' (^) ' [A 2 PC + A 2 PV ] (29) 

Here uncertainties connected with the weak interaction vertex are emphasized and 
our choice for the corresponding most uncertain parameters is shown explicitly. We 
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recall once again that the estimate ( |2"8"D depends on the validity of the factorization 
approximation, which in particular, allows to neglect the photon emission from final 
state quarks. 

In order to compare our resluts for the WA mechanism to the short distance 
(SD) penguin contribution, we cast the latter into the form of eq. (p9|). Using the 
results of |5[ [J we obtain 

| A PC (B- -> p- 7 ) SD | = | A PV (B- -> p- 7 ) SD | 
-0.12GeV-f^«VH) . (30) 



Thus, our result in eq. (|22|) amounts to about a 10% correction to the SD ampli- 
tude for B~ — > p~7 with uncertainty at the level of 50% from the values of the 
CKM matrix elements. This correction is comparable to the present accuracy of the 
hadronic matrix elements in refs. |5], ||; with more refined techniques, however, the 
WA contributions should presumably not be neglected in predicting the branching 
ratio of B~ — > p~7 and the CP- violating asymmetries ||. In comparison, the phe- 
nomenological estimates of the long distance penguin contributions 0, O] indicate 
that these are quite small and probably negligible. The correction of WA mechanism 
in B° — > p° 7 turns out to be much smaller, at the level of 1%, due to relative factor 
a 2 /ai and due to relative smallness of the amplitudes ( p3| ) as compared with fl2"2"|). 

Unlike in radiative S-decays, the WA effects are crucial in the corresponding D- 
decays. The short distance penguin contributions are completely negligible, while 
the long distance terms may even lead to observable rates [11], |32|, |33|| : 

BR{D+ - p+ 7 ) = ^(V cd V: dai ff p ml {^^)\^pc + A 2 PV ]r D - (31) 



Using the values in eqs. (|24D and (p2q), we obtain 

BR(D + -> p+ 7 ) ~ 2.7 • 10~ 6 , BR(D° -> p° 7 ) ~ 3.1 • 10~ 6 (a 2 /0.5) 2 . (32) 



With the relations (27) we predict the branching ratio for D° — > K*° 7 to be 



1.5 • 10" 4 (a 2 /0.5) 2 and for D s -> p+ 7 to be 2.8 • 10~ 5 . 



We can apply the present technology to the decays B — > /z/^ 7 [34]; we simply 



replace the p current with the leptonic one. With the notation of ref. 3"5|] , we obtain 

^ = 247r2 i 2m 4 f 1 ^(-1 - P 2 ) V{A%c(p 2 ) + A 2 PV (P 2 )} (33) 

for the ratio between the purely leptonic decay modes with and without photon. A 
problem arises because we now need the form factor for all values of p 2 and not just 
at the p mass, and the sum rule method becomes unreliable for large p 2 values. We 
can overcome the difficulty by assuming a pole behaviour (with B* pole for Ape and 
axial B pole for Apy) and normalizing it with our numbers obtained from the sum 
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rules (|2"DD and (PT|) at p 2 = lbGeV 2 . The phase space suppression at large values of 
the lepton pair invariant mass makes the result insensitive to the accuracy of this 
extrapolation. We then obtain ~ 19.0 which is in agreement with the earlier 
work Q and underlines the importance of the radiative leptonic decays. 

Finally, we should mention that the same correlation function ((5|) may be used 
to derive the sum rule for the B*B^f and D*D^i couplings. The derivation follows 
essentially the same steps as in where the light-cone sum rules for the B*Btt 
and D*Dtt couplings have been obtained. The detailed analysis is beyond the scope 
of the present investigation and will be presented in a separate paper p7| . We just 
mention here that our preliminary results on the photonic D* width combined with 
its pionic widths obtained in |26j agree well with experimental data. 



7. In this paper we have calculated a part of the long-distance contributions to 
exclusive radiative rare B and D decays. We have shown that the QCD sum rules 
with the light-cone photon wave function allow to calculate them rather reliably. In 
the case of the decays B —>■ p7 the long-distance contributions considered are small, 
especially in the neutral mode , which allows to neglect them in the determination 
of V t d through that decay. We predict branching ratios for the corresponding D 
decays where the long-distance weak annihilation dominates. There are other non- 
perturbative effects which can be estimated in a similar fashion; work on this is in 
progress. 
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(a) (b) 

Figure 1: Weak annihilation mechanism for the decay B~ — > p~~{. Two additional diagrams 
with photon emission from the final state quarks are not shown. 




(d) (e) (f) 

Figure 2: Diagrams corresponding to OPE of the correlation function in eq.(6): (a), (d) 
perturbative contributions, (b) quark, (c) quark-gluon condensate contributions; (e) emission 
of photon at large distances parametrized by photon wave function; (f) gluon correction to 
the diagram (e). Solid lines denote quarks, dashed lines -gluons, wavy line -photon, arrows - 
external currents. Crosses denote vacuum quark and gluon condensates, ovals in the diagrams 
(e) and (f ) denote the photon light-cone wave function. 



